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We have investigated the effects of quantum fluctuations of quasiparticles on the operation of 
superconducting radio-frequency single-electron transistors (RF-SETs) for large values of the quasi- 
particle cotunneling parameter a = 8Ej/E c , where Ej and E c are the Josephson and charging 
energies. We find that for a > 1, subgap RF-SET operation is still feasible despite quantum fluctu- 
ations that renormalize the SET charging energy and wash out quasiparticle tunneling thresholds. 
Surprisingly, such RF-SETs show linearity and signal-to-noise ratio superior to those obtained when 
quantum fluctuations are weak, while still demonstrating excellent charge sensitivity. 
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The radio-frequency single electron transistor (RF- 
SET) is a highly sensitive, fast electrometer, and has been 
suggested as a potentially quantum-limited linear am- 
plifier suitable for measurements of individual electronic 
charges 0, S 0] • Recent investigations have addressed 
use of the RF-SET as an electrometer^ 0, @ , a readout 
device for charge based qubits 0, H, @, and a sensor for 
real-time electron counting experiments |l0|. Linearity 
is a fundamental assumption of theoretical discussions 
of the quantum limits of amplifiers QUO]- Nonetheless, 
there has been no detailed investigation of the range of 
linear response for the RF-SET. 

Most theoretical studies of RF-SET performance focus 
on normal metal SETs, either in the sequential tunnel- 
ing H2I Il3| or cotunneling regimes , while most ex- 
periments are performed using a superconducting SET 

(SSET) HEIIEE3- Trans P° rt in thc SSET can be di- 
vided into two regimes, depending on the relative sizes of 
the bias voltage and superconducting gap A: above- 
gap (eVdc > 4A), dominated by Coulomb blockade of 
quasiparticles, and subgap (eVdc < 4A), dominated by 
combinations of quasiparticle and resonant Cooper pair 
tunneling known as Josephson-quasiparticle (JQP) cycles 
[I4I Hif . While the best charge sensitivities are found for 
above-gap operation 0, the SSET backaction — the rate 
at which it dephases a measured system — is largest there 
0, E| . Recent work has focused on subgap operation 
for which backaction is significantly reduced, and shot 
noise is non-Poissonian [j, @, Uij . Theoretical studies of 
quantum fluctuations in the SSET have been limited to 
above-gap cotunneling of quasiparticles 01 ■ In this Let- 
ter we find that linearity and subgap quantum charge 
fluctuations in superconducting RF-SETs are intimately 
related: as quantum fluctuations strengthen, linearity 
and signal-to-noise ratio (SNR) improve, while charge 
sensitivity remains excellent. 

Our SSETs consist of a small island connected to 
macroscopic leads via two Al/AlO^/Al tunnel junctions 
Jl(2) with normal state resistances i?i(2) as illustrated 
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FIG. 1: (a) Schematic diagram of the SET illustrating RF op- 
eration. A voltage Din consisting of dc and RF biases Vdc and 
v r f is incident on a tank circuit consisting of an inductor L, a 
capacitor C p , and the SET, with tunnel junction resistances 
and capacitances Ri(2) an ^ ^1(2) ■ A small charge oscillation 
qo cos u) m t modulates the reflection coefficient F of the tank 
circuit and therefore the reflected voltage v r . (b) Electron 
micrograph of S2 (taken after all measurements). Gates Gl 
and G2 were used vary the SET offset charge, (c) Power spec- 
trum of v r for qo = 0.063e rms and cu m /2n = 100 kHz. The 
measured sideband power and noise floor were used to find 
the charge sensitivity and SNR of the RF-SET. 



in Fig.^a). We have fabricated and characterized three 
samples, SI, S2, and S3 with total resistance R n — 
R\ + R2 of 58, 38 and 24kf2; an electron micrograph of 
S2 is shown in Fig. ^b) . The samples were mounted on 
the mixing chamber of a dilution refrigerator at its base 
temperature of 20 mK. A Nb chip inductor L w 120 nH 
together with the parasitic capacitance C p « 0.2 pF of 
the SET contacts constituted a tank circuit with reso- 
nant frequency /lc ~ 1 GHz and quality factor Q w 16. 
We measured the samples' current-voltage (I-V) char- 
acteristics in an asymmetric voltage-biased configuration 
[Fig. Ola)] by varying the dc bias voltage Vd c in the ab- 
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TABLE I: Sample parameters. Resistances are in kQ, energies 
in fieV, and areas in 10~ 3 fim 2 . 
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sence of an RF excitation. Details of RF operation are 
similar to those discussed elsewhere jj, ||j . The SET off- 
set charge Qo+qo cosuj m t consisted of a dc component Qo 
that set the overall working point and an ac component 
of amplitude qo that modulated the reflected voltage v r . 
Power spectra of v r [Fig. [Jc)] were used to determine 
the charge sensitivity 5q and SNR. 

In Fig. [3 we show representative I-V characteristics of 
the samples in the superconducting state, measured for 
different Qo, with qo = 0. For SI, we observe clear above- 
gap (Vdc 800 /xV) current modulation corresponding to 
Coulomb blockade of quasiparticle tunneling [Fig.^a)]. 
The main sub-gap features corresponding to the JQP 
01 cycles are sharp and clearly distinguished. As 
illustrated in Fig. the simplest JQP cycle consists of 
resonant tunneling of a Cooper pair through one junction 
and dissipative tunneling of two quasiparticles through 
the other, transporting two electrons through the SET. 
The cycle can occur only when the transition — ► 1 
(1 — > 0) is allowed, i. e., for eVdc > E c + 2A where 
E c = e 2 /2Cs is the charging energy of the SET and 
Cs = Ci + C2 + 2C ff its total capacitance. While the JQP 
cycle is forbidden at lower bias, at Qo/e = n g w i and 
eKic = 2E C Cooper pair tunneling is resonant at both 
junctions and the double JQP (DJQP) cycle becomes 
possible. The fact that sequential tunneling cannot oc- 
cur via either cycle for 2E C < eVdc ^ E c + 2 A is reflected 
in SI by a sharp drop in current at Vdc ~ 630 fiV just 
below the JQP feature. 

As R n decreases, so does current modulation for 
eVdc > 4A, consistent with suppression of the Coulomb 
blockade by quasiparticle cotunneling ^3 : the modula- 
tion is reduced for S2, and nearly absent for S3 [Fig.^b) 
and (c)]. In contrast, features corresponding to the JQP 
cycles still exist but become progressively less sharp. 
Since these cycles involve both Cooper pair and quasi- 
particle tunneling, we hypothesize that subgap quantum 
fluctuations of quasiparticles are strong, while quantum 
fluctuations of Cooper pairs remain weak. Since to the 
best of our knowledge no theoretical description of sub- 
gap quantum charge fluctuations in the SSET exists, we 
provide simple arguments supporting our hypothesis. 

We first compare with known results for above-gap 
transport. Following Ref. Jj| we define a parameter 
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a = J;ff(i?r 1 + i ?2 X ) = TT characterizing the strength 
of quantum fluctuations for quasiparticles, assuming 
^1(2) = and using the Ambegaokar-Baratoff rela- 
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FIG. 2: I-V characteristics for (a) SI (b) S2 and (c) S3 
(note scale change), were chosen for Qo showing the DJQP 
process (red), the JQP process (blue), and an intermediate 
value of Qo (green). The arrows and vertical hash mark show 
the peak-to-peak RF amplitude 2Qv T { and dc bias Vdc for 
optimal RF-SET operation. Inset: variation in the measured 
charging energy E c relative to the bare charging energy E® 
for SI (solid triangle), S2 (circle) and S3 (square). Error bars 
indicate uncertainty in E®. Solid line: theoretical prediction. 



tion for the Josephson coupling energy Ej = — -jp where 
Rk = ■ Quantum fluctuations are negligible for a <C 1 . 
Determining E c from the location of the DJQP peak and 
Ej from the total junction resistance we calculate a as 
in Table |U None of our samples satisfies a€l, although 
for SI (a = 0.78) some above-gap Coulomb modulation 
survives. The progressively weakening modulation for S2 
(a = 1 .08) and S3 (a = 2.65), is consistent with previous 
results |l7j . 

Cotunneling as described in Ref. [13 occurs only 
for Vdc > 4A/e: it results in two quasiparticle ex- 
citations and transfers a single electron through the 
SET. Other virtual processes, however, remain impor- 
tant for Vdc < 4A/e. For normal SETs, E c can be 
renormalized by quantum charge fluctuations: e. g., 
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(a) JQP (b) DJQP (c) VJQP 

FIG. 3: Various JQP cycles. Here J2(l) is on the left (right) 
and Vdc > 0. Solid (empty) circles indicate quasielectrons 
(quasiholes) created during a cycle, (a) JQP cycle. Beginning 
in the state n = (n = 1), where n is the number of excess 
electrons on the SET, the transition — > 1 (1 — > 0) is allowed, 
bringing Josephson tunneling through J 1(2) into resonance. 
Cooper pair tunneling 1 — 1 (0 <4> 2) is interrupted by 
quasiparticle tunneling through the opposite junction — 1 — » 

(2 — > 1), completing the cycle, (b) DJQP cycle. When 
Josephson tunneling is simultaneously resonant through both 
Jl and J2, transport occurs via the sequence -O- 2, 2 — ► 1, 

1 — 1, —1 — » 0. (c) Proposed VJQP cycle. If the transition 
—* 1 (1 — » 0) is forbidden, it may still occur virtually. The 
remaining JQP transitions are allowed for relevant Vdc. 



near n g — 0, the effective charging energy E c « 
E®(1 — 4g) where g — Rk/^Ru is the dimension- 
less parallel conductance of the tunnel junctions and 
E® the bare charging energy; similar renormalization 
occurs in the superconducting state || 0|. Calcu- 
lating the first-order energy shift due to transitions 
n — > nil, we find the renormalized charging en- 
ergy E* = E° c (l + g^T[^(l + 2n B )] + T[^(l - 2n g )}) 

where T(x) = J Q °° K^{u)e~ xu du and Ki(u) is a Bessel 
function. 

Using the expression for E^, we find empirically that 
E® = 254 /xeV gives the measured E c for SI. We 
measure the total geometric junction area j4 t ot for the 
samples with an estimated accuracy of ±20%, obtain- 
ing the values in Table [I] Setting E® = e 2 /2C^ where 
C° = C? + Cl + 2C g and using 2C g « 80 aF, we obtain 
Ci+C® — 195aF as the total unrenormalized junction ca- 
pacitance for SI. Scaling this result according to A tot we 
find C% E° c and finally E s c for S2 and S3 [Tabled; agree- 
ment is excellent given the uncertainties in At Q t- The 
inset to Fig. [21 shows the relative difference between E c 
and E® scaled by 1/g. The results agree with theory to 
within our experimental accuracy, providing strong evi- 
dence that subgap quantum fluctuations of quasiparticles 
occur in our samples. 

Virtual quasiparticle tunneling may also play a role in 
subgap transport, as suggested by the softening of the 
JQP cycle cutoff in S2 and S3. To illustrate such effects 
more clearly we show a plot of the /(Vd c ,%) surface for 
S2 in Fig. ||(a). The JQP resonances along the <^ 2 
and 1 — 1 lines and the DJQP peak at their intersec- 
tion are clearly visible, but there is no sharp cutoff of 
the JQP process below the 1 — ► (0 — > 1) thresholds. 



For comparison, in Fig. 0(b) we show a simulation of the 
current in S2 based on sequential tunneling |l9| at an ele- 
vated temperature and including photon-assisted tunnel- 
ing due to an electromagnetic environment. Despite the 
extreme conditions the quasiparticle tunneling thresholds 
are clearly visible, and the SSET current drops nearly to 
zero between the JQP and DJQP features. The absence 
of quasiparticle thresholds in Fig. 0(a) calls for an expla- 
nation outside of the sequential tunneling picture. 

A candidate process that could allow transport along 
the Cooper pair resonance lines between the JQP and 
DJQP features is illustrated schematically in Fig. EJc). 
If below threshold the transition 1 — > (0 — ► 1) occurs 
virtually, the transitions <^> 2 and 2 — ^ 1 ( — 1 1 and 
1 — > 0) are allowed, completing what we call the virtual 
JQP (VJQP) cycle. Two quasiparticle excitations are 
created, but two electrons are transferred through the 
SET, so that the process should be allowed for eV > 2A. 
The energy barrier Ej, for 1 — > (0 — * 1) vanishes at 
threshold and climbs to Eb « E c + 2 A at the DJQP peak. 
The process can be neglected if the allowed quasiparti- 
cle tunneling rate T qp is small compared to the inverse 
dwell time of the virtual quasiparticle: T qp <C Eb/K. Us- 
ing T qp = 4A/e 2 i?„, this becomes R n > which 
is violated for a range of voltages between the DJQP 
and JQP features. A detailed theoretical analysis is re- 
quired to determine the contribution of the VJQP cycle 
to transport. 

In contrast to the quasiparticle thresholds, features as- 
sociated with Cooper pair tunneling are visible in both 
the data and the simulation, suggesting that the num- 
ber of Cooper pairs is well defined. For the JQP pro- 
cess at resonance, the Coop er pair tunneling rate is 
r cp w Ej/hT qp = f "X- Hot- Demanding that energy 
broadening due to Cooper pair tunneling be small com- 
pared to the typical energy barrier AE C for virtual tun- 
neling gives 2Sr cp /4£ , c = jf-^f- <§C 1, which is easily 
satisfied even for S3. For S2 and S3, then, quantum fluc- 
tuations are significant for quasiparticles but small for 
Cooper pairs. 

We now turn to RF operation. Optimal operating con- 
ditions were selected as follows: a small charge oscillation 
qo ~ 0.006e rms was applied and the SNR determined 
from the power spectrum of v r as in Fig.^c) . Subgap op- 
eration (all samples) and above-gap operation (SI) were 
optimized over dc bias Vd c , rf bias v r { and offset charge 
Qo- We measured SNR versus input amplitude qo for 
each optimization and determined the charge sensitivity 
Sq using Sq — , q ° io~ SNR / 20 where the resolution band- 

° VBW _ 

width BW = 1 kHz and SNR is in dB 0. 

The optimized biases for SI and S3 are indicated in 
Fig. [5] and the results of the Sq and SNR measurements 
in Fig. □ For SI the best Sq = 9 x 10~ 6 e/y/Hz was 
found for Vdc = 860 /iV, consistent with previous re- 
sults [||. Linearity, however, was poor: as go increases, 
the measured SNR rapidly becomes sublinear, and Sq 
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FIG. 4: False color images of I(Vd c ,n g ) for (a) S2 at 
T — 20 mK (b) a simulation at T = 200 mK assuming an 
electromagnetic environment with impedance i? C nv = 50 f2 
and temperature T cnv = IK. Cooper pair resonance lines 
<4> 2 ( — 1 <^> 1) and quasiparticle tunneling thresholds 1 — *• 
(0 — ► 1) are indicated by the dashed and solid lines. 



worsens [Fig. Et a )]- Since Sq apparently does not sat- 
urate even for go = 4.5 x 10 _3 e rms it is unclear how 
small go must be to achieve linear response. For sub- 
gap operation (Vd c = 600 fiV) of SI [Fig. ED 3 )], we 
find Sq ps 1.3 x 10~ 5 e/VHz, with SNR nearly linear to 
go ^ O.Ole rms. Since Sq appears close to saturation at 
go = 3.1 x 10~ 3 e rms, we may have approached linear 
response. 

For S3 the best operating point occurred at Vd c = 
440 /zV [Fig. Etc)], between the DJQP and JPQ features 
with Sq w 1.2 x 10 ~ 5 e/VSz, better than that for sub- 
gap operation of SI. Moreover, linearity was vastly im- 
proved: the SNR remains linear and Sq nearly flat to 
go = 0.038e rms indicating that we have achieved linear 
response in this sample. For S2 (data not shown) the 
best Sq ~ 1.2 x 10~ 5 e/VHz also occurred subgap, and 
the SNR was linear to go w 0.02erms. 

We can now make some general statements about the 
effects of quantum fluctuations on RF-SET operation. 
For samples with smaller a such as SI, transport is fairly 
well described by the sequential tunneling picture: I-V 
characteristics are sharp and vary strongly with Qq giv- 
ing rise to excellent charge sensitivity. The same sharp- 
ness, however, prevents good linearity, since a large go 
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FIG. 5: Charge sensitivity Sq and SNR (linear scale) versus 
go in e rms for (a) SI, above gap, (b) SI, subgap, and (c) 
S3, subgap. Charge sensitivity (solid red symbols) is plotted 
on the left axis and SNR (open blue symbols) on the right. 
For reference, the SNR for linear response is plotted as the 
dashed lines for Sq measured at the smallest go. 



necessarily moves the SET far from optimal operation. 
For samples with larger a such as S3, quantum fluctu- 
ations cause at least two important effects. First, the 
subgap features are smoothed and broadened, improving 
linearity: e. <?., in S3 there is no "dead spot" between 
the DJQP and JQP features for which the SSET current 
is roughly independent of Qo- Second, renormalization 
of E c moves the DJQP feature to lower bias, so that the 
optimal rf amplitude of about (2A — E c )/e increases with 
a. Finally, the smaller R n simplifies impedance matching 
between the RF-SET and the 50 coaxial line. 

In conclusion, we have investigated the influence of 
quantum charge fluctuations on the charge sensitivity 
and SNR of RF-SETs. We find that RF-SETs with 
a > 1-2 (strong quantum fluctuations) show both good 
linearity and good charge sensitivity. In contrast, RF- 
SETs with a < 1 (weak quantum fluctuations) show 
poor linearity and only modestly better charge sensitiv- 
ity. These findings assume particular importance given 
interest in the RF-SET as a potentially quantum-limited 
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linear amplifier. We have achieved linear response only 
for subgap operation in samples with a > 1 for which 
quantum fluctuations of quasiparticles are substantial. 
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